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Motivated by recent experimental interest in tunneling into heavy electron materials, we present 
a theory for electron tunneling into a Kondo lattice. The passage of an electron into a Kondo lattice 
is accompanied by a simultaneous spin flip of the localized moments via cotunneling mechanism. 
We compute the tunneling current with the large- N mean field theory. In the absence of disorder, 
differential tunneling conductance exhibits two peaks separated by the hybridization gap. Disorder 
effects lead to the smearing of the gap resulting in a Fano lineshape. 

PACS numbers: 71.27.+a, 75.20.Hr, 74.50.+r 



Major developments in scanning tunneling electron 
spectroscopy (STEM) over the last decade,_particularly 
as a probe of cuprate superconductors [E [1, H HJ] , suggest 
that this tool will find increasing utility as an atomic- 
scale probe of many-body phenomena in new classes of 
materials. One area of particular promise lies in the ap- 
plication of STEM to heavy fermion materials. 

Heavy fermion compounds contain a dense lattice of lo- 
calized magnetic moments interacting with a sea of con- 
duction electrons to form a "Kondo lattice" [H, 0] . These 
materials exhibit a diversity of many body behaviors, in- 
cluding anisotropic superconductivity, Kondo insulating 
behavior and quantum criticality. Motivated by recent 
tunneling experiments on /-electron materials [3, B 
in this paper we develop a theory for tunneling into a 
coherent Kondo lattice. 

How do electrons tunnel into a Kondo lattice, where 
the main degrees of freedom are local moments? Since di- 
rect tunneling into localized magnetic orbitals is blocked 
by Coulomb interactions, the naive expectation is that 
the electrons can only tunnel into the surrounding con- 
duction sea. In 1960s Anderson and Appelbaum (lol . 
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111 . [12| recognized that magnetic ions actively participate 
in the tunneling process via a "cotunneling mechanism" 
[l3L 14 1 in which the passage of a tip electron into the 



conduction sea occurs cooperatively with a spin-flip of 
localized moments. The manifestation of cotunneling in 
the tunneling conductance of quantum dots and mag- 
netic atoms adsorbed on surfaces is well established ex- 
perimentally 13, 3, 15, lij] . Here we examine the effect 
of these processes on tunneling into a coherent band of 
excitations of a Kondo lattice, deriving a new expression 
for the tunneling current into a Kondo lattice in terms 
of the Green's function of composite co-tunneling oper- 
ators. Using the large-A^ approximation, we show how 
cotunneling processes open a direct tunneling channel be- 
tween the tip and the composite quasiparticle states of 
the Kondo lattice. Once coherence develops, cotunnel- 
ing and direct tunneling processes interfere, giving rise 
to distinctive two peak structures in tunneling spectra. 

We begin by writing down the Kondo lattice Hamil- 
tonian in the presence of a tunneling probe, which takes 
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FIG. 1: (Color online) Electron tunneling into a heavy- 
fermion material involves two parallel processes: direct tun- 
neling with amplitude t c into the conduction sea and cotun- 
neling with amplitude t / into a composite combination of the 
conduction electron and local magnetic /-moments. These 
composite states are expected to develop coherence below the 
Kondo temperature Tk- Inset shows a typical differential 
conductance curve observed for tunneling into a single Kondo 
ion. 



the form H = H KL + H, 



tip 



Ht, where 



HKL = e k c L c k<r + ^/C?) ' (cja^CjT?) (1) 

k,cr j 

is the unperturbed Kondo lattice Hamiltonian, Cj a — 
c k(Te ik ' Rj creates a conduction electron and of(j) 
is the spin operator of a localized /-electron at site j, re- 
spectively. The term H t i P = X)k<r £ kPk a Pi<.a describes the 
electrons in the tip. The crucial new feature of this model 
lies in the composite character of the tunneling Hamil- 
tonian. When the tip lies in the vicinity of site 0, the 
tunneling Hamiltonian is given by 



where 



Ht = PoalpOa + H.c, 
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contains a direct tunneling term of amplitude t c and a 
"cotunneling term" of amplitude tf. From the equations 
of motion, the tunneling current operator is 



as 118 



eN r = 



POa 



H.c. 
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where N c — J^ko- °La c ^ ^ s the number operator of the 
conduction electrons. From the form of I and Ht, we 
see that the passage of an electron from the tip into the 
lattice is accompanied by a spin- flip of a local moment. 
In this way, one particle states in the tip are coupled to 
the composite fields which define heavy electron quasi- 
particles. 

The Hamiltonian ((2]) is a Kondo lattice generaliza- 
tion of the Anderson- Appelbaum tunneling Hamiltonian 



111 [12| , first introduced to explain zero-bias anomalies 
associated with tunneling between two metallic leads via 
a single localized moment. Similar models have subse- 
quently been used to describe tunneling through a quan- 
tum dot [14j . The cotunneling component of Ht can 
be understood as a result of mixing between states in 
the tunneling tip and the localized orbitals of the Kondo 
lattice. This process distorts the symmetry of the Wan- 
nier states that hybridize with the localized moments, 
partially delocalizing them into the tip. A derivation 
of the cotunneling terms can be done by carrying out a 
Schrieffer- Wolff transformation on the Anderson model 
describing the lattice and the tip [12, 0, In the 

Anderson model, the localized /-electrons hybridize with 
the conduction electrons. When a tip is introduced above 
site of the lattice, tunneling between the /-state and 
the probe electrons modifies the hybridization according 
to H h — > (Vcl a + tfpl a )f 0cr + H.c, where tf is the am- 
plitude to tunnel directly from an /-state to the probe, 
so that the tip modifies the orbital hybridizing with the 
/-state: 



CQa -> C 0cr + yPOcr- 



(5) 



After a Schrieffer- Wolff transformation is carried out, 
which reduces the Anderson model to a Kondo model, 
this same replacement must be made to the Kondo inter- 
action at site in the unperturbed Kondo Lattice model. 
To leading (linear) order in tf/V, the result of this pro- 
cedure is the quoted result in ([3]), where tf = Jtf/V. 

Next, we compute the tunneling current. One of the 
questions that immediately arises, is whether the differ- 
ential conductance can be analyzed in a conventional way 
when cotunneling terms are present. We now show that 
even though / contains a composite operator, the weak- 
ness of the tunneling matrix elements still permits us to 
expand the current to leading order in the tunneling ma- 
trix elements, thereby rewriting it in terms of the full 
many-body Green's functions of the bulk. To carry out 
this procedure, we write steady-state tunneling current 



I(eV) ee (I) = ^Re / ^G* ( W ), 
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where G^(w) is the Keldysh Green function [lj| between 
the tip electron and 4>0a ■ Expanding the current to lead- 
ing order in the tunneling matrix elements, we obtain 
O = G*G$ + G%G§, where R, A, K denote the 
retarded, advanced and Keldysh Green's functions of the 
tip and the Kondo lattice. Since the tip and the lead 
are in thermal equilibrium, their Keldysh Green's func- 
tions can be re- written in terms of retarded and advanced 
Green's functions, using the fluctuation dissipation re- 
lations [lj] Gf(u) = -2inp p (Lj + eV)h{uj + eV) and 
G^(lo) = — 2itt p^(uj)h{ijj) . Here h(iv) = 1 — 2/(w), where 
f(w) is the Fermi distribution function, while pu P {u)) 
is the local density of states of the tip; p^{uj) is the 
"cotunneling" density of states of the sample given by 
p^{to) = ^ImG,/, (a; — i5), and G^(w) is the retarded 
Green's function of the ip field, usually obtained through 
analytic continuation of the Matsubara imaginary time 
propagator Gy, = -{Tip 0o ,(r)i/jl a (0)). 

Using these relations, the current ([6]) can be re-written 



I(eV) =2s£/dw 



(/(«- 



- eV)pip(w) 
eV) - /H). 
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In this way, the tunneling current into a Kondo lattice 
probes the spectral function of the composite operator. 

To illustrate the tunneling into the Kondo lattice, we 
now solve for the tunneling behavior in the large- N limit 
[IE HE HlHi of the Kondo lattice, where N = 2j + l 
is the spin degeneracy of the localized /-state. In this ap- 
proach, the spin operator is represented as a bilinear of 
pseudo-fermions [2J], Sf(j) = fja.S a pfjp, where S a p are 
generators of the SU (N) symmetry group. The mean- 
field theory provides a representation of the composite 

fcrmion [a a p ■ Sf(J) \ Cjp in ([3]) as a single fermionic op- 
erator 



(fap ■ Sf{j)) c, , 



V 



(8) 



where the amplitude ^ = —{fjpCjp}- In this way, the 
large- N mean field theory captures the formation of a 
composite /-electron, an essential element of the Kondo 
effect. In terms of pseudo-fermions, we can re- write single 
particle operator in ([3]) as 



Yjot — t c Cj a -(- tffj a7 



(9) 



where the complex amplitude for tunneling into the com- 
posite fermion state is tf — j-tf. 

The requirement that the number of pseudo-fermions 
at any given site should be equal to N/2 introduces a 
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constraint A, to be determined self-consistently together 
with the hybridization amplitude V (see e.g. 
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The resulting mean-field Hamiltonian can then be di- 
agonalized by means of the Bogoliubov transformation 
Ckcr = v k a ka + u k b kcr , and / kCT = u k a kcr - v k b ka , 
where u k and v k are the Kondo lattice coherence fac- 

?k 

The Hamiltonian |T|) in 



tors given by u 2 . = [R k + (e k - A)]/2i? k , vt = 1 - " 



with i? k = v /(e k -A) 2 +4V 2 . 
the mean-field approximation then becomes H^J 1 ^ = 

Dka( w k a ka^kQ+^k& kQ &k Q ), where = (e k +A±i? k )/2 
is the quasiparticle dispersion in the newly developed 
heavy Fermi liquid. The mean-field tunneling Hamilto- 
nian then becomes 



4 m/) =E4 UcCja + tff: 



+ H.c. 



(10) 
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Although our mean field Hamiltonian has the form of 
the Anderson lattice model with U — 0, the states on 
which it operates have an underlying composite struc- 
ture, formed when local spins hybridize with conduc- 
tion electrons. Thus, the Hamiltonian (flO|) provides a 
mean-field description of the tunneling into the conduc- 
tion band together with the cotunneling processes involv- 
ing local moments. 

It is instructive to contrast the tunneling conductance 
expected in a Kondo lattice with that of a single Kondo 
impurity. Using the tunneling Hamiltonian (|10[) , we com- 
pute G^(tu). In the case of a single Kondo impurity, we 
obtain 



(t c iirpV + i 



- A 



(11) 



where p is the density of states of the conduction elec- 
trons, A = irpV 2 ~ Tk is the width of the Kondo reso- 
nance. The differential conductance jjk = g(eV) is 
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9im P (eV) = N —r—t\ Ptip P 



\q- 



l + e' ; 



=(eV-A)/A 



(12) 

where N is the spin degeneracy and q — A(eV)/B(eV) 
is the ratio of two complex tunneling amplitudes, where 
A(eV) = if+t c VP( eV 1 _ ) describes the cotunneling into 
the atomic orbital and B(eV) = trVTr5(eV — e k ) describes 
direct tunneling into the metal [26|. Here e' = (eV — 
A)/ A, pup is the density of states at the Fermi level of 
electrons in the tip. For a broad flat band, A = if, 
B = t c Virp and q = tf/(t c Vnp). 

Now we turn to the case of the Kondo lattice. Within 
the large- N mean field theory, we obtain 



G« L (u;) = Nj2 



(tc + tfz^:) 2 



£k 



V 2 
jj — A 



(13) 



where e k is the dispersion of the conduction band. We 
obtain the following expression for the differential tun- 



neling conductance, 



l(eV) = N 



2ire 2 



s=±,k s k 

(14) 

where i? sk = (w sk — A)/A. The prefactor of the delta 



function has a characteristic Fano functional form [271 . 
|28| . This form introduces an asymmetry in the resulting 
voltage dependence of the tunneling conductance g(eV). 
The momentum summation in G^ L (ui) (fT3|) and g(eV) 
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FIG. 2: (Color online) Differential tunneling conductance 
g(V) for a single Kondo impurity case (blue line) given by 
(|12[) . a Kondo lattice (green line) given by Q5) . A typical 
Fano shape in the single Kondo impurity case gets replaced 
with a double-peaked resonance line in the Kondo lattice case. 
The dashed lines illustrate the effect of disorder, which de- 
stroys the coherence, closing the gap in the density of states 
curve. Here t f /t c = .2, q = 4.9, \/T K = .3, while D = W0T K . 



(|14p can be carried out analytically assuming a constant 
conduction electron density of states p, to give 

'27re 2 \ „ 1. 



g(eV) = N 



where 



tip tip p -1mG" h (eV-iS), (15) 



qA 



(u) = (H — r) In 

1 — A 



L UJ — A 

uj — \ 



Do 



2D/t 2 c 



iv — A 
(16) 

Here —D\ and D2 are the lower and the upper conduction 
band edges respectively, and 2D = D\ + D2 is the band- 
width. The differential tunneling conductance predicted 
by this formula has two well-pronounced peaks at eV ~ A 
separated by a narrow hybridization gap A g ~ 2V 2 /D in 
the single particle spectrum, as shown in Fig. 2. 

In practice, experimental tunneling results will be 
modified by the effects of disorder [20(. A phenomeno- 
logical quasiparticle elastic relaxation rate T may be in- 
troduced into the theory by replacing u> — > to — iT in 
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([13|) . The results of this procedure are shown in Fig. 
2. As we see, disorder removes the sharp peak structure 
in the tunneling conductance g(eV) (f!5|) . The resulting 
lineshape of the tunneling conductance dl / dV (eV) is an 
asymmetric smooth curve. 

The current work can be extended in a number of in- 
teresting directions. One important aspect, is to examine 
the effects of cotunneling on the fluctuations in the den- 
sity of states probed in Fourier transform STM exper- 
iments. In one-band systems, the Fourier transform of 
these fluctuations is phase sensitive to quasiparticle scat- 
tering 0, 0| , and is expected to be an important probe 
of both the quasiparticle dispersion and the phase of the 
cotunneling matrix elements. 

A particularly fascinating aspect of cotunneling is its 
likely interplay with various forms of heavy fermion or- 
der, such as heavy fermion superconductivity. Unlike in 
conventional tunneling, the quasiparticle matrix elements 
of the composite operators associated with cotunneling 
are expected to be sensitive to the nature of the heavy 
electron ground-state. For example, recent work has pro- 
posed that heavy electron superconductivity may involve 
composite pairing between local moments and electron 
pairs . A key feature of composite pairing is the pres- 
ence of two conduction screening channels Ti and T2, so 
that now the tunneling will be described by the field 
© of the form 



uf a +v sgn(a)/l c 



IpOa = t c C 0a + ^2 



i=l 



i/r, (oWs • 5/(0)) cr i/3 



(17) 



where v describes hybridization in the particle-particle 
channel. In this way, we see that the cotunneling term 
in ip may develop both particle and hole components, 
resulting in Andreev reflection even in the limit of weak 
tunneling. 

In conclusion, wc have studied electron tunneling into 
a Kondo lattice of localized moments, bringing out the 
importance of cotunneling as a primary mechanism of 
tunneling into the heavy electron fluid. We have ex- 
pressed the conductance in terms of a spectral function 
of a cotunneling composite operator, illustrating the re- 
sult by a calculation carried out in the large- N limit. 
Our results predict that in a clean system the differential 
tunneling conductance will display two peaks separated 
by the hybridization gap. Addition of disorder leads to 
the smearing of the gap and produces a Fano-like smooth 
asymmetric lineshape. 
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